Secondary Structures in Long Compact Polymers

Compact polymers are self-avoiding random walks which visit every site on a lattice. This polymer model is used widely for studying statistical problems inspired by protein folding. One difficulty with using compact polymers to perform numerical calculations is generating a sufficiently large number of randomly sampled configurations. We present a Monte-Carlo algorithm which uniformly samples compact polymer configurations in an efficient manner allowing investigations of chains much longer than previously studied. Chain configurations generated by the algorithm are used to compute statistics of secondary structures in compact polymers. We determine the fraction of monomers participating in secondary structures, and show that it is self averaging in the long chain limit and strictly less than one. Comparison with results for lattice models of open polymer chains shows that compact chains are significantly more likely to form secondary structure.Comment: 14 pages, 14 figure

Dynamical fluctuations in biochemical reactions and cycles

We develop theory for the dynamics and fluctuations in some cyclic and linear biochemical reactions. We use the approach of maximum caliber, which computes the ensemble of paths taken by the system, given a few experimental observables. This approach may be useful for interpreting single-molecule or few-particle experiments on molecular motors, enzyme reactions, ion-channels, and phosphorylation-driven biological clocks. We consider cycles where all biochemical states are observable. Our method shows how: (1) the noise in cycles increases with cycle size and decreases with the driving force that spins the cycle and (2) provides a recipe for estimating small-number features, such as probability of backward spin in small cycles, from experimental data. The back-spin probability diminishes exponentially with the deviation from equilibrium. We believe this method may also be useful for other few-particle nonequilibrium biochemical reaction systems

Unbiased sampling of globular lattice proteins in three dimensions

We present a Monte Carlo method that allows efficient and unbiased sampling of Hamiltonian walks on a cubic lattice. Such walks are self-avoiding and visit each lattice site exactly once. They are often used as simple models of globular proteins, upon adding suitable local interactions. Our algorithm can easily be equipped with such interactions, but we study here mainly the flexible homopolymer case where each conformation is generated with uniform probability. We argue that the algorithm is ergodic and has dynamical exponent z=0. We then use it to study polymers of size up to 64^3 = 262144 monomers. Results are presented for the effective interaction between end points, and the interaction with the boundaries of the system

Sequence Heterogeneity Accelerates Protein Search for Targets on DNA

The process of protein search for specific binding sites on DNA is fundamentally important since it marks the beginning of all major biological processes. We present a theoretical investigation that probes the role of DNA sequence symmetry, heterogeneity and chemical composition in the protein search dynamics. Using a discrete-state stochastic approach with a first-passage events analysis, which takes into account the most relevant physical-chemical processes, a full analytical description of the search dynamics is obtained. It is found that, contrary to existing views, the protein search is generally faster on DNA with more heterogeneous sequences. In addition, the search dynamics might be affected by the chemical composition near the target site. The physical origins of these phenomena are discussed. Our results suggest that biological processes might be effectively regulated by modifying chemical composition, symmetry and heterogeneity of a genome.Comment: 10 pages, 5 figure

A lattice model of hydrophobic interactions

Hydrogen bonding is modeled in terms of virtual exchange of protons between water molecules. A simple lattice model is analyzed, using ideas and techniques from the theory of correlated electrons in metals. Reasonable parameters reproduce observed magnitudes and temperature dependence of the hydrophobic interaction between substitutional impurities and water within this lattice.Comment: 7 pages, 3 figures. To appear in Europhysics Letter

Evolution of the potential-energy surface of amorphous silicon

The link between the energy surface of bulk systems and their dynamical properties is generally difficult to establish. Using the activation-relaxation technique (ART nouveau), we follow the change in the barrier distribution of a model of amorphous silicon as a function of the degree of relaxation. We find that while the barrier-height distribution, calculated from the initial minimum, is a unique function that depends only on the level of distribution, the reverse-barrier height distribution, calculated from the final state, is independent of the relaxation, following a different function. Moreover, the resulting gained or released energy distribution is a simple convolution of these two distributions indicating that the activation and relaxation parts of a the elementary relaxation mechanism are completely independent. This characterized energy landscape can be used to explain nano-calorimetry measurements.Comment: 5 pages, 4 figure

Nonuniversal power law scaling in the probability distribution of scientific citations

We develop a model for the distribution of scientific citations. The model involves a dual mechanism: in the direct mechanism, the author of a new paper finds an old paper A and cites it. In the indirect mechanism, the author of a new paper finds an old paper A only via the reference list of a newer intermediary paper B, which has previously cited A. By comparison to citation databases, we find that papers having few citations are cited mainly by the direct mechanism. Papers already having many citations ('classics') are cited mainly by the indirect mechanism. The indirect mechanism gives a power-law tail. The 'tipping point' at which a paper becomes a classic is about 21 citations for papers published in the Institute for Scientific Information (ISI) Web of Science database in 1981, 29 for Physical Review D papers published from 1975-1994, and 39 for all publications from a list of high h-index chemists assembled in 2007. The power-law exponent is not universal. Individuals who are highly cited have a systematically smaller exponent than individuals who are less cited.Comment: 7 pages, 3 figures, 2 table

A Bell-Evans-Polanyi principle for molecular dynamics trajectories and its implications for global optimization

The Bell-Evans-Polanyi principle that is valid for a chemical reaction that proceeds along the reaction coordinate over the transition state is extended to molecular dynamics trajectories that in general do not cross the dividing surface between the initial and the final local minima at the exact transition state. Our molecular dynamics Bell-Evans-Polanyi principle states that low energy molecular dynamics trajectories are more likely to lead into the basin of attraction of a low energy local minimum than high energy trajectories. In the context of global optimization schemes based on molecular dynamics our molecular dynamics Bell-Evans-Polanyi principle implies that using low energy trajectories one needs to visit a smaller number of distinguishable local minima before finding the global minimum than when using high energy trajectories

Model for Folding and Aggregation in RNA Secondary Structures

We study the statistical mechanics of RNA secondary structures designed to have an attraction between two different types of structures as a model system for heteropolymer aggregation. The competition between the branching entropy of the secondary structure and the energy gained by pairing drives the RNA to undergo a `temperature independent' second order phase transition from a molten to an aggregated phase'. The aggregated phase thus obtained has a macroscopically large number of contacts between different RNAs. The partition function scaling exponent for this phase is \theta ~ 1/2 and the crossover exponent of the phase transition is \nu ~ 5/3. The relevance of these calculations to the aggregation of biological molecules is discussed.Comment: Revtex, 4 pages; 3 Figures; Final published versio